function L = LEs(f, df, param, x0, Ndiscard, Ncompute)
% LEs 计算一维离散映射的李雅普诺夫指数（单个）
%
%   L = LEs(f, df, param, x0, Ndiscard, Ncompute)
%
%   输入参数：
%       f        - 映射函数句柄，形如 f(x, param)
%       df       - 映射的导数函数句柄，形如 df(x, param)
%       param    - 控制参数(可为标量，也可为向量)
%       x0       - 初始值
%       Ndiscard - 舍弃的前期迭代步数
%       Ncompute - 用于计算李雅普诺夫指数的迭代步数
%
%   输出参数：
%       L        - 对应的一维映射的李雅普诺夫指数近似值

    % 1) 先舍弃前 Ndiscard 步
    x = x0;
    for i = 1:Ndiscard
        x = f(x, param);
    end

    % 2) 正式计算李雅普诺夫指数
    sumLog = 0;
    for i = 1:Ncompute
        x = f(x, param);
        d = df(x, param);
        sumLog = sumLog + log(abs(d) + 1e-15);
    end

    % 3) 平均
    L = sumLog / Ncompute;
end
